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Three-dimensional Dirac and Weyl semimetals exhibit a disorder-induced quantum phase transi¬ 
tion between a semimetallic phase at weak disorder and a diffusive-metallic phase at strong disorder. 
Despite considerable effort, both numerically and analytically, the critical exponents v and 2 of this 
phase transition are not known precisely. Here we report a numerical calculation of the critical 
exponent v = 1.47 ± 0.03 using a minimal single-Weyl node model and a finite-size scaling analysis 
of conductance. Our high-precision numerical value for v is incompatible with previous numerical 
studies on tight-binding models and with one- and two-loop calculations in an e-expansion scheme. 

We further obtain z = 1.49 ± 0.02 from the scaling of the conductivity with chemical potential. 


Introduction. Materials with an electronic band- 
structure dispersing linearly from a Fermi point are 
among the driving themes in contemporary condensed 
matter physics w- After the experimental verifi¬ 
cation of such a Dirac-type bandstructure in single¬ 
layer graphene [JJ, the focus has now turned to three- 
dimensional materials. The compounds NaaBi and 
CdaAs 2 have been confirmed as “Dirac semimetals” 0- 
0. In materials that break either time- or space inver¬ 
sion symmetry, the twofold band degeneracy of Dirac 
semimetals is lifted and the resulting phase is termed 
Weyl semimetal. The non-centrosymmetric compounds 
TaAs and NbAs have recently proven experimentally to 
harbor such Weyl nodes in their bandstructures dDKGp. 
Similar bandstructures have been achieved in a photonic 
crystal realization in Ref. Ha- 

Theoretical work accompanied and, in part, preceded 
the recent experiments. Beyond the single particle pic¬ 
ture, Coulomb interactions were argued to be marginally 
irrelevant in the renormalization group (RG) sense due 
to the vanishing density of states at the Fermi point 
mm- Quenched disorder, however, inevitably present 
in realistic materials, is a much more subtle issue. Dat¬ 
ing back to work from the 1980s on mi, the presence 
of a disorder-induced quantum phase transition is by 
now firmly established analytically mmm and nu¬ 
merically EMI- In the weak-disorder phase, the ran¬ 
dom potential is irrelevant in an RG sense. Thus, for 
large system sizes and low temperatures, a weakly dis¬ 
ordered system qualitatively behaves as a clean system 
with renormalized Fermi velocity. This leads to a num¬ 
ber of experimentally important predictions for weak dis¬ 
order, such as quadratically vanishing density of states 
mm or pseudoballistic charge transport m at the 
nodal point. In contrast, for strong disorder one finds a 
metallic phase with finite density of states at the Fermi 
energy and diffusive transport characteristics [21, 27 . 

Signatures of the disorder-induced quantum critical¬ 
ity are expected in almost any experimentally relevant 
observable, from heat capacity to transport properties. 
Standard scaling theory |2S| predicts power-law depen¬ 
dences on disorder, chemical potential, or temperature 
in the vicinity of the critical point mm- The only 
input to this variety of predicted power laws is a pair of 


critical exponents characteristic of the universality class. 
Denoting the dimensionless disorder strength and chem¬ 
ical potential by K and p, respectively, close to the criti¬ 
cal point K = K c , p, = 0 the correlation length exponent 
v and the dynamical critical exponent z govern the rela¬ 
tion between reduced disorder strength k = \K — K c \/K c 
and the emerging correlation length ^ as £ oc k~ v , and 
the relation between emergent energy- and length scales 
as e oc (~ z . Although the critical point is located at zero 
chemical potential, predictions of scaling theory persist 
for small finite doping. 

To date, the best analytical estimates for v and z for a 
single Weyl- or Dirac node follow from a Wilsonian mo¬ 
mentum shell RG calculation in an e-expansion scheme 
around critical dimension two. The results of the one- 
loop calculation by Goswami and Chakravarty are v = 1 
and z = 1.5 [3]. The accuracy of the one-loop exponents 
was challenged by a calculation of two-loop diagrams 
by Roy and Das Sarma ESI, who found v = 1.14 and 
z = 1.31. On the other hand, there are instances where 
the e expansion strategy is known to fail completely, the 
Anderson metal-insulator transition in three dimensions 
being a well known example m — although the present 
transition is of a different type as it connects two non¬ 
insulating phases EH:- Numerical results for the critical 
exponents obtained from tight-binding models harboring 
multiple Weyl- or Dirac nodes E2H2B] are in reasonable 
agreement with the one-loop results above, albeit with 
large uncertainties in v , z and K c . 

Motivated by the lack of a firm theoretical prediction 
and in view of potential experiments, we performed a 
numerical calculation of the critical exponents in a sin¬ 
gle Weyl node using state of the art finite-size scaling 
for quantum transport properties. Our results, which 
we report in detail below, have significantly reduced un¬ 
certainties in comparison to the previously known es¬ 
timates. Whereas our result for the dynamical critical 
exponent, z = 1.49±0.02, is consistent with the previous 
numerical calculations and with the one-loop calculation 
(but not with the two-loop calculation!), our value for 
the correlation length exponent, v = 1.47 ± 0.03, devi¬ 
ates rather significantly. 

Minimal model and numerical method. The minimal 
model for the disorder induced quantum criticality is a 


2 


single Weyl node with potential disorder, 

H = hvcr ■ k + p + U(r ), (1) 

where v is the Fermi velocity, <x denotes the vector of 
Pauli matrices, and k measures the Bloch wavevector rel¬ 
ative to the nodal point. We connect the Weyl semimetal 
to two ideal leads, both modeled as Weyl nodes with 
p taken to infinity and without the random potential 
U. The Weyl semimetal has dimension 0 < x < L 
and 0 < y, z < W in transport and transverse direc¬ 
tions, respectively. To quantize transverse momenta 
k y<z we apply periodic or antiperiodic boundary condi¬ 
tions (PBC/APBC). An ultraviolet cutoff A restricts the 
magnitude of transverse wavevector \k V}Z \ < A and 1/A 
sets the microscopic length scale. The random potential 
U (r) is assumed to have zero mean and Gaussian white 
noise fluctuations 

(U(r)U(r')) = X(hv) 2 6(r-v'), (2) 

with I\ the dimensionless disorder strength. The chem¬ 
ical potential p has to vanish to reach the critical point, 
however, we will work with finite p below to assess the 
dynamical critical exponent z. 

We study the signatures of disorder-induced quantum 
criticality of Eq. ([l]) in a quantum transport framework 
at zero temperature, employing the numerical scattering 
matrix method of Ref. [21], which is based on related 
studies of disordered Dirac fermions in two dimensions 
on Eg. The conductance can be computed from the 
scattering matrix’s transmission block t using the Lan- 
dauer formula G = trtU and is measured in units of e 2 /h 
throughout. 

Correlation length exponent v: Finite-size scaling of 
the conductance for p = 0. The standard method to 
assess the correlation length exponent v is finite-size 
scaling m- To perform such an analysis, one needs to 
identify a dimensionless observable that assumes differ¬ 
ent values on the two sides of the (bulk) phase tran¬ 
sition. In Ref. m we showed numerically that the 
conductance G fulfills these requirements. For large as¬ 
pect ratio r = W/L 1 and in the thermodynamic 
limit L — > oo, the conductance takes the values Gk= o = 
r 2 (ln2/27r) ~ O.llr 2 in the pseudoballistic phase at dis¬ 
order strength K = o m\ and Gk^too = <Jr 2 L —> oo in 
the diffusive phase for I\ > K c (with bulk conductivity 
a) m- In the vicinity of the critical point K = K c , when 
the system dimensions L, W are larger than all internal 
length scales other than the emerging correlation length 
(, G assumes a scaling form G = G(L/(/,W/Q. Using 
C oc k~ v and fixing the aspect ratio r = W/L, we arrive 
at G = G r {L 1 ' v k), with universal correlation length ex¬ 
ponent v and a scaling function G r that depends on r 
and the boundary conditions. 

Numerically, the conductance G is found to vary con¬ 
siderably for different disorder realizations, however with 
the restriction G > Gk=o for every disorder realiza¬ 
tion. The histogram of SG = G — Gk =o is shown in 
Fig. [I] for the specific choice L = 27r/(Ar) x 11, PBC. 


For weak disorder, K < 5, we find that the distribu¬ 
tion of SG can be well fitted by a log-normal distribu¬ 
tion p^jsiiSG) = e _ i ln5G ^ MLN i 2 / 2<T LN/ ( 7 iAr ^G' v / 27 r ! with 
parameters /xln and ctln, see Fig. |T] A feature not cap¬ 
tured by this fit is the tail of large but rare conduc¬ 
tances, which are possibly related to rare region effects 
m- As G is not self-averaging (analogous to Ander¬ 
son localization), we choose the median m of the dis¬ 
tribution as a scaling quantity, i. e., we search a scaling 
function m = m r {Lfl v F). The standard error of the 
median is calculated using the asymptotic variance for¬ 
mula = l/4p(m) 2 N, where N is the total number of 
disorder realizations and the unknown exact probability 
density p is approximated by a smooth interpolation of 
the measured histogram. For further research, it would 
be desirable to understand the occurrence of the empir¬ 
ical log-normal conductance distribution. 



g—g k=0 


Figure 1: (color online) Probability distribution of the differ¬ 
ence SG = G — Gk=o of the conductance with and without 
disorder potential. The aspect ratio r = W/L = 5, periodic 
boundary conditions were applied in the transverse direction, 
and the sample length L = (2n/A) x 11/r. The dimensionless 
disorder strength is indicated in the figure. The number of 
disorder realizations is 4000. Dots indicate the value of the 
median m, which is used as the scaling variable. Dashed lines 
show fits to a log-normal probability distribution. The data 
for K > 4.25 is offset in vertical direction for clarity. 


We compute m for a range of disorder strengths K 
and lengths L, for aspect ratios r = 5 and 7, and we 
also varied the boundary condition between periodic and 
anti-periodic. The data for r = 5, PBC, is shown in Fig. 
[2] the other data sets can be found in the supplemental 
information, Ref. [ 35 ] . At criticality, where £ diverges, 
m r is independent of L and the data traces in Fig. [2]cross 
in one point. In the supplemental information [35] . we 
show the details of a least squares fit for m r ( L 1 ' I 'k ) for 
small k to a polynomial of fourth order in Lfl v k (solid 







3 



K 

Figure 2: (color online) Scaling plot for the logarithm of the 
median m of the distribution of 8G = G — Gk=o for a dis¬ 
ordered Weyl node with W/L = 5 and periodic transverse 
boundary conditions. For better visibility we plot log(m) — K 
vs. K on the vertical axis. The solid curves show the results 
of a least-squares fit to a scaling form m = m r {L 1 ' v k), ex¬ 
panded in a fourth order polynomial (for details, see 13511. 
The estimates for the most important fit parameters and their 
standard error are K c = 5.024±0.004, log(m c ) = 0.12±0.004, 
v = 1.47 ± 0.02. The quality of the fit is X 2 /N = 0.93. The 
sample lengths are indicated in the figure. The gray vertical 
line indicates the position of the estimated critical disorder 
strength. The inset shows a scaling collapse of the data in 
the main panel using the estimated values of K c and v. 

lines). An excellent and stable fit was achieved even 
without including any irrelevant scaling variable that we 
took in leading order as L v with y < 0. Taking into 
account the fitting results of all other parameter sets in 
a standard procedure (see [35]) we find v = 1.47 ± 0.03. 
The conductance data for smaller aspect ratios r < 3 
(data not shown) reveals a large irrelevant contribution 
to the scaling function that hindered a successful fit in 
terms of a simple low order polynomial. 

In Ref. [2T] it was argued that the Fano factor F (the 
ratio of shot-noise power and conductance) is an alter¬ 
native quantity to distinguish the pseudoballistic from 
the diffusive phase. In the pseudoballistic phase one 
has F(K < K c ) ~ 0.574 while in the diffusive phase 
F(K > K c ) = 1/3. In the Supplemental Material we 
show that our result for v is consistent with the value ob¬ 
tained from a finite-size scaling analysis using the Fano 
factor (v = 1.40 ± 0.05). This method however suffers 
from an inferior quality of the data set, both in terms of 
error bars of the individual data points as well as in the 
range of system sizes available. 


Dynamical critical exponent: Scaling of critical bulk 
conductivity. We now turn to the dynamical critical ex¬ 
ponent z that connects the emergent length scale £ and 


the corresponding energy scale e in the vicinity of the 
fixed point. In our transport geometry, a natural choice 
of a quantity that has a scaling involving the dynamical 
exponent z is the bulk conductivity a, which is also of 
immediate experimental relevance. 

To connect the dynamical critical exponent to the con¬ 
ductivity, we again start with a scaling form around crit¬ 
icality m■ Since the unit of a in three dimensions is 
inverse length, we find a (fc, y) = £ _1 /(/z/e) with an un¬ 
known dimensionless scaling function /. We define a new 
scaling function / ( x ) = x 1 ^ f (a; -1 ) in terms of which 
a(k,y) = y 1 ^ z f(k zv /y) . At K = K c , the ‘critical’ 
conductivity tr c thus scales as 

cr c (m) k 1/z - (3) 

The scaling form (J3| is valid with small corrections 
within an extended quantum critical region [20] for fi¬ 
nite k when the argument of / is sufficiently small, i.e. 
k -C k* (y) oc y}'^ uz \ This allows us to numerically 
compute an estimate of 0 in spite of the fact that the 
value of K c is known only within error bars. 

We compute G(L) for fixed large W, PBC, a range 
of y and for K = 5.0, which is within the K c confidence 
interval m■ We perform a disorder average over at least 
ten disorder realizations. Since transport in a Weyl node 
at finite y is diffusive, we expect G = <jW 2 /L, which 
is confirmed in the simulation. Finite size effects are 
irrelevant once W, L are larger than the characteristic 
//-induced length scale oc y~ x / z . We show tr c vs. y for 
I\ = K c in Fig. [3] (dots) and indeed observe a power 
law (solid line) for y < Hvh. with inverse exponent z = 
1.49±0.02. For larger chemical potentials, y comparable 
to the band edge Hv A, the scaling breaks down and from 
Drude transport theory we expect a crossover to a oc y 2 , 
proportional to the density of states (dashed line). 


Discussion. We numerically studied the disorder- 
induced quantum phase transition in three-dimensional 
Dirac materials in terms of a minimal model, a sin¬ 
gle Weyl node with potential disorder. In contrast to 
the well-known Anderson metal-insulator transition, this 
disorder-induced phase transition for a single Weyl node 
is between two non-insulating phases. In addition to the 
correlation-length exponent iq it features a non-trivial 
dynamical critical exponent z, which has no counterpart 
in the standard metal-insulator transition. 

Our high-precision results for the exponents v and z 
not only allow for a variety of quantitative predictions of 
experimentally observable power laws around criticality 
— such as the density-of-states exponent /3 = v(3 —z) for 
k > k c m but also improve on previously reported 
theoretical predictions. Our result v = 1.47±0.03 differs 
significantly from analytical results obtained from a one- 
or two-loop £-expansion RG calculation (y = 1, 1.14, re¬ 
spectively nun!])- The failure of the e-expansion cal¬ 
culation is reminiscent of the situation for the Anderson 
localization in three dimensions |3(Jj , where v = 1.375 in 
the symplectic class [3B]. Our estimate for the dynamical 
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Figure 3: Scaling of the conductivity a(K,y) for (near) 
critical disorder strengths K = 5.0 (dots) with chemical po¬ 
tential /r. The solid line is a power-law fit a oc to the 
data points, with 2 = 1.49 ± 0.02. The dashed line indi¬ 
cates a it oc (i 2 power law as expected from Drude trans¬ 
port theory when scaling breaks down. The sample width 
W = (2-7r/A) x 29 and periodic boundary conditions were 
applied. For chemical potentials below the values shown in 
the figure a bulk conductivity could not be reliably obtained 
from the calculated conductance data. 


critical exponent is z = 1.49 ± 0.02, in agreement with 
the one-loop RG calculation (z = 1.5), but not with the 
two-loop prediction (z = 1.31). 

In principle, the model of Eq. ([!]), which has white- 
noise disorder and a sharp momentum cutoff, could be 
modified to include a more faithful representation of the 
microscopic disorder, albeit at an increased numerical 
cost. For example, Ref. m employed a finite disor¬ 
der correlation length that sets the microscopic length 
scale; the mode cutoff then can safely be taken to infinity. 
However, the difference between these two models is ir¬ 
relevant in the RG sense and thus both models are in the 
same universality class. To see this, recall that a finite 
disorder correlation length is equivalent to a higher-order 
momentum dependence of the disorder-induced interac¬ 
tion vertex in the replicated disorder averaged action and 
is thus irrelevant m ■ On the other hand, the numerical 
value of K c is non-universal and, thus, sensitive to the 
disorder model. In this context we note that a model 
with sharp momentum cut-off has also been used in Ref. 
m, where it was found to give the same results as a 
model with finite disorder correlation length. Moreover, 
in a realistic band structure the linear form of Eq. ([l]) is 
only an approximation. Quadratic corrections, however, 
are RG irrelevant and thus will not change the critical 
exponents m- 

Realistic Weyl and Dirac semimetals. Realistic Weyl 
or Dirac semimetals have multiple Weyl nodes [35], ei¬ 
ther separated in momentum space or distinguished by 
their transformation properties under point group sym¬ 
metries. The same applies to numerical studies based 
on tight-binding models [22H25], which confirmed the 


presence of a disorder-induced phase transition on the 
basis of density-of-states calculations. With multiple 
Weyl nodes, disorder might not only cause intra- but 
also inter-node scattering of quasiparticles. The latter 
process is not captured by our minimal model. Symme¬ 
tries in more realistic models with multiple Weyl nodes 
may also be different from the minimal model: While our 
minimal model has an effective time-reversal symmetry 
mapping the single Weyl node onto itself, in realistic 
models, time-reversal or inversion symmetries can relate 
different nodes or be absent. Although the precise na¬ 
ture of disorder potentials in realistic three-dimensional 
Dirac materials is yet to be determined, there are plau¬ 
sible scenarios in which intra-node scattering dominates 
over inter-node scattering, a priori justifying the use of 
our minimal model. For example, in Weyl semimetals 
the ratio of scattering rates is controlled by the smooth¬ 
ness of the disorder potential and the separation of Weyl 
nodes in momentum space m- 

The applicability of the minimal model in the pres¬ 
ence of sizeable inter-node scattering, i. e., the ques¬ 
tion whether or not the presence of some amount of 
inter-node scattering changes the universality class of 
the disorder-induced semimetal-metal transition, is an 
issue that has not been conclusively settled [10] . Inter¬ 
node scattering is omitted in field-theoretical approaches 
PH US]- Empirical evidence that inter-node scattering 
does not affect the universality of the transition comes 
from Ref. [21] . which found a remarkable universality 
for three different disorder types in a tight-binding Dirac 
semimetal model, albeit with large error bars on the crit¬ 
ical exponents. 

If the assertion of a single universality class insensitive 
to inter- or intra-cone scattering (and the related sym¬ 
metry differences) is correct, the observed critical expo¬ 
nents should match with those obtained in tight-binding 
models. However, such studies [221125] yield a value for 
v inconsistent with our result, for a typical example see 
Ref. |22j which finds v = 0.92 ± 0.13 at K > K c and 
v = 0.81 ± 0.21 at K < K c . The value of the dynam¬ 
ical critical exponent z = 1.5 ± 0.1 from Ref. m is 
in agreement with our result. Assuming that the type 
of scattering is indeed immaterial for critical exponents, 
we attribute the large difference with the tight-binding 
model exponent v to difficulties in accurately estimating 
the critical disorder strength from density-of-states data. 
The uncertainty of K c translates to a large uncertainty 
in the critical exponent v. In contrast, our very precise 
estimate of K c was possible using the finite-size scal¬ 
ing method where K c can be obtained from the unique 
crossing of the data in Fig. [2] 

In the supplemental information [35] , we exemplify 
this interpretation by revisiting the density of states data 
from Ref. [22] (cf. Fig. 3a) obtained at zero energy for 
a range of disorder values around the critical disorder 
strength. Using the critical exponent /3 = 2.22 calcu¬ 
lated with our estimates for v and z we are able to pro¬ 
duce an excellent fit for the data points in the vicinity 
of the critical disorder strength, though we find a much 
smaller critical disorder strength than asserted in Ref. 
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[22] ■ Since the microscopic model in Ref. [22] and in this 
work are different, the values of the non-universal critical 
disorder strengths cannot be compared. Uncertainty in 
K c does not cause comparable problems when determin¬ 
ing the critical exponent z, because the large size of the 
critical region in the chemical potential-disorder param¬ 
eter plane renders the extraction of z much less sensitive 
to the uncertainty in K c . This is consistent with the mu¬ 
tual agreement between our estimate for z and the value 
in Ref. [22] , 
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Supplemental Material 


Details of the finite-size scaling analysis. We here provide details of the finite-size scaling procedure, following 
Refs. pmFiT| . In addition to the data set presented in the main text — aspect ratio r = 5 and periodic boundary 
conditions (PBC) —, we have obtained conductance distributions for antiperiodic boundary conditions (APBC) with 
r = 5 and for aspect ratio r = 7, PBC. 

The sample width W = L ■ r is set to be W = (2n/A )(M — 1/2) for APBC and W = (2n/A)M for PBC, with M 
a positive integer. The transverse wavenumbers are k y>z = (2n/W)n VtZ , with n VtZ = —M, —M + 1,..., M for PBC 
and n ViZ = —M ± 1/2, — M + 3/2,..., M — 1/2 for APBC. A summary of all data sets used in this work is given in 
Table [D 

For each data set, the median m(K , M) of the conductance distribution is determined. We perform a least-squares 
fit to a polynomial of the form 


m(K, M) = ao(l + b 01 L y + ... + b 0qo L yyo ) + ai ■ (L^ v k) • (1 + b n L y + ... + b lqi L vqi ) (4) 

+ • • • + a p ■ (L^kY • (1 + b n L y + ...b lqp L yy *) 


for medians obtained at the same value of the aspect ratio r and the same boundary conditions. Data points (i.e., 
medians of conductance distributions) and fits are shown in Fig. [2] of the main text for r = 5 and PBC, and in Fig. 
[4] for r = 5, APBC, and and r = 7, PBC. The following algorithm for the fitting procedure is used: The order of the 
polynomials in Eq. Q is increased by adding a new fit parameter m or if (i) the merit function \ 2 /N £ [0, oo] 
(N is the number of data points) for the resulting fit is lowered by more than 2% compared to the previous fit and 
(ii) the error of any fitting parameter (as calculated from error propagation theory) does not exceed the parameter’s 
estimate in magnitude. Initial values for each fitting procedure are chosen randomly and the parameter estimates for 
the best fit out of a few hundred fitting trials is reported in Table [T] along with the error estimates and the value of 
X 2 /N. A fit is acceptable if y 2 /A r < 1, another measure is the ‘goodness of the fit’ Q £ [0,1] where Q = 1 indicates 
a perfect fit (for definitions of Q and y 2 /iV see, for example, Ref. [43]). 

Ideally, fitting parameters should not strongly depend on the number of different values of M within a data set. 
We successfully checked the stability of the fitting results by repeating the fitting procedure above for reduced data 
sets (deleting data points of the largest or smallest M in the data set r = 5, PBC), as indicated in Table [TJ Finally, 
the estimate for v is calculated as an average of the best fit estimates for v for each data set whereas the total error 
bars are unions of error bars from each single data set (’practical-error-bar procedure’, see Ref. [H]V 


r = W/L 

B.C. 

M 

N 

x 2 /n 

G 

V 

K c 

log{m c ) = a 0 

ai • 10 1 

a 2 ■ 10 2 

—a 3 ■ 10 2 

—a 4 • 10 3 

5 

PBC 

11,16,21,26,31 

70 

0.93 

0.4 

1.47 ±0.02 

5.024 ±0.004 

0.12 ±0.004 

5.0 ±0.1 

3.0 ±0.3 

1.3 ± 0.1 

1.8 ±0.9 

5 

PBC 

16,21,26,31 

56 

0.62 

0.94 

1.46 ±0.02 

5.007 ±0.007 

0.102 ±0.006 

4.9 ± 0.2 

3.3 ±0.4 

1.2 ± 0.1 

2.6 ±1.2 

5 

PBC 

11,16,21,26 

56 

0.84 

0.55 

1.48 ±0.02 

5.036 ±0.006 

0.128 ±0.004 

5.0 ±0.1 

3.2 ±0.4 

1.4 ± 0.1 

2.3 ±1.3 

5 

APBC 

14,19,24,29 

56 

0.97 

0.29 

1.47 ±0.02 

5.031 ±0.005 

0.169 ±0.004 

4.9 ± 0.1 

2.8 ±0.3 

1.1 ±0.1 

1.7 ±0.8 

7 

PBC 

19,26,33 

42 

0.64 

0.87 

1.47 ±0.03 

4.983 ±0.009 

0.782 ±0.007 

4.8 ± 0.2 

2.1 ±0.2 

1.5 ±0.2 

- 


Table I: Details of the finite-size scaling procedure. The left part of the table specifies the data sets subject to a least squares 
fit with model (|4]| while the right part gives the fitting results. Numbers with ± are error bars (one standard deviation). The 
range of disorder strength for all data sets is K = 4.25, 4.375, ..., 5.875. 
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K K 

Figure 4: (color online) Scaling plots for the logarithm of the median m of the distribution of the difference 5G = G — Gk= o 
of the two-terminal conductance with and without disorder potential, for a Weyl node with aspect ratio r = W/L = 5 and 
antiperiodic boundarity conditions (left) and r = 7 and periodic boundary conditions (right). For clarity, we plot log(m) — K 
vs. K. The solid curves show the results of a least-squares fit to a scaling form m = m r expanded in a fourth (r = 5) 

and third (r = 7) order polynomial, respectively. The estimates for the most important fit parameters are given in Table |I] 
The gray vertical line indicates the position of the estimated critical disorder strength. 


Comparison with finite-size scaling for Fano factor. As discussed in the main text, besides the conductance, 
also the Fano factor can be expected to be a suitable observable for a finite-size scaling analysis. A scaling plot is 
shown in Fig. [5] and the details of the analysis (done as above for the conductance data) are reported in Table [TT) 
Although the number of disorder realizations is comparable to the corresponding conductance data in Fig. [4] (left), 
for the Fano factor error bars are much larger. Moreover, while for conductance scaling data traces for system sizes 
M = 14,19,24,29 all cross in a single point, the Fano factor data for M = 14 does not cross with the traces of the 
larger system sizes, indicating that shot noise around criticality is controlled by larger emergent length scales than 
conductance. For the remaining system sizes, the analysis yields v = 1.40 ± 0.05. Given the intrinsic difficulties for 
the Fano factor data discussed above we consider the error bar overlap with the conductance result v = 1.47 ± 0.03 
as a confirmation for consistency of the two finite-size scaling methods. 


W/L = r 

B.C. 

M 

N 

x 2 /n 

s 



log(m c ) = a 0 

ai • 10 1 

a 2 • 10 2 

a 3 ■ 10 3 

a 4 ■ 10 3 

5 

APBC 

19,24,29 

42 

0.72 

0.69 

1.40 ±0.05 

4.994 ±0.015 

-3.4 ±0.01 

4.3 ±0.5 

5.9 ± 1.5 

-6 ±2 

-5 ±3 


Table II: Details of the finite-size scaling procedure for the Fano factor. The left part of the table specifies the data sets subject 
to a least squares fit with model @ while the right part gives the fitting results. Numbers with ± are error bars (one standard 
deviation). The range of disorder strength is K = 4.25, 4.375, ..., 5.875. 
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Figure 5: (color online) Scaling plots for the logarithm of the median m of the distribution of the difference SF = Fk =o — F 
of the two-terminal Fano factor without and with disorder potential, for a Weyl node with aspect ratio r = W/L = 5 and 
antiperiodic boundarity conditions. For clarity, we plot log(m) — K vs. K. The solid curves show the results of a least-squares 
fit to a scaling form m = m r expanded in a fourth order polynomial. The data set for M = 14 was not included in 

the analysis. The estimates for the most important fit parameters are given in Table |TT] The gray vertical line indicates the 
position of the estimated critical disorder strength. 


Comparison with density-of-states scaling for tight-binding model. We revisit the results of a recent density- 
of-states simulation in a disordered Dirac semimetal from Ref. |22j . The study is based on a large four-band 
tight-binding model tuned at the topological phase transition between a strong and weak topological insulator. If 
inter-node processes can be neglected, around criticality the density of states at zero energy should increase as 
p(e = 0) oc (K — K C )P with /3 = (3 — z)v. Using our estimate (3 = 2.2, we successfully fit the data from from Ref. 
[22| . (cf. Fig. 3a) in Fig. [b] (solid line), except for the three data points with largest disorder strength. In contrast, 
the emphasis in the interpretation of Ref. m was laid on data points for larger K , excluding the immediate vicinity 
of the critical point at K = K c . This leads to a larger estimate for K c and a smaller estimate for (3. 



Figure 6: (color online) Bulk density of states of a disordered tight-binding model of a Dirac semimetal at zero energy as a 
function of disorder strength around criticality. The data (dots) is adapted from Ref. |22| . Fig. 3a. The solid line is a fit to 
the scaling form p[t = 0) oc (K — K c ) /3 with (3 = (3 — z)v fixed to 2.22 from our finite-size scaling analysis taking into account 
data points with K < 7 only. The vertical gray line shows the corresponding estimate K c = 6.03. The dashed line is a fit of 
the eight data points with the largest I\ with setting K c = 6.4 [22), the estimate for the exponent is (3 = 1.16. 






